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I.   INTRODUCTION 

In  the  high  energy  laser  field,  the  electric  discharge  convection  laser 
(EDCL)  offers  good  efficiency  as  well  as  simplicity  of  handling  relative  to 
gasdynamic  and  chemical  lasers.   Progress  towards  higher  laser  outputs  from 
EDCLs,  however,  has  been  impeded  by  the  collapse  of  the  glow  discharge  at  the 
increased  currents  and  pressures  desired.   This  glow  collapse  (usually  re- 
ferred to  as  arcing  or  constriction)  represents  the  practical  upper  limit  of 

the  electrical  power  that  may  be  pumped  into  the  necessary  non-equilibrium 
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processes.   It  has  been  shown    that  aerodynamic  stabilization  can  be  used 

very  successfully  in  overcoming  the  glow-collapse  limits.   It  should  be 
pointed  out  here  that  the  laser  application  is  but  one  application  of  the 
relatively  new  interdisciplinary  field  of  cold  plasma  chemistry. 

Aerodynamic  stabilization  is  based  on  two  separate  effects  usually  re- 
ferred to  as  convection  (laminar)  and  mixing.   The  purely  convective  effect 
is  said  to  decrease  the  residence  time  of  the  particles  through  the  discharge 
and  thereby  to  counteract  unstable  mechanism  that  lead  to  the  glow  collapse  . 
The  mixing  effect,  which  may  be  due  to  vortex  flows  or  to  turbulence  from 
boundary  layers  or  grids,  has  been  observed  to  have  powerful  stabilizing 
characteristics.   Convection  scales  up  with  flow  velocity  whereas  mixing 
depends  on  the  intensity,  spectral  distribution,  and  location  of  the  eddies 
or  vortices. 


II.   SPECIFICS 

At  the  Technion,  the  author  had  the  opportunity  to  participate  in 

experimental  work  dealing  with  vortex-flow  stabilization  '   and  with  super- 
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sonic  flow  stabilization  '  .   These  techniques  complement  the  work  that  has 

been  underway  at  the  Postgraduate  School  on  turbulence  stabilization  of 

discharges.   Some  results  of  the  Technion  work  have  already  appeared  (see  Refs, 

6-11)  and  more  will  be  forthcoming. 

A  significant  portion  of  the  author's  time  was  dedicated  to  a  study  of 
theoretical  aspects  of  discharge  stabilization.   The  author  was  privileged  to 
colaborate  with  Dr.  Y.  Khait  in  establishing  a  perspective  regarding  some  of 
the  author's  previous  results  with  turbulence  stabilization.   This  effort 
resulted  in  the  writing  of  a  paper  titled  "Influence  of  Turbulence  on  Diffuse 
Electrical  Gas  Discharges  Under  Moderate  Pressures,"  which  is  being  submitted 
for  publication.   Other  theoretical  work  involved  a  close  look  at  the 
Schottky  solution  for  the  positive  column  in  a  glow  discharge  with  the  aim 
at  establishing  its  use  for  a  stability  formulation.   This  effort  was  fruit- 
ful on  two  intermediate  aspects  of  the  problem  and  these  are  included  here 
as  Appendixes  A  and  B. 

A  proposal  was  written  to  the  United  States-Israel  3inational  Science 
Foundation   which  extensively  outlines  the  progress  of  the  entire  effort. 
If  awarded,  the  cooperation  of  the  Technion  and  the  Postgraduate  School 
would  be  formalized  (however,  the  budget  would  provide  neither  salaries  for 
principal  investigators  nor  equipment  for  the  USA  experiments) . 


III.   DISCUSSION 

The  classical  solution  to  the  ambipolar  diffusion  of  charges  in  the 
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positive  column   *   of  a  glow  discharge  offers  no  hints  as  to  how  striations 

and  constriction,  among  other  important  observations,  manifest  themselves  in 

the  discharge.   As  pressure  and/or  current  are  increased,  instabilities  grow 

in  the  discharge  and  it  is  important  to  have  a  sufficiently  complete  model 

of  the  physics  which  reflects  these  effects.   In  other  words,  before  flow 

effects  such  as  (laminar)  convection  and  turbulence  are  introduced  it  is 

necessary  to  re-examine  some  assumptions  inherent  in  the  Schottky  solution. 

A  partial  list  of  phenomena  not  reflected  in  the  Schottky  solution  may 

be  given  as  follows: 

1)  effect  of  sheaths  at  the  walls 

2)  effect  of  electrodes  and  intermediate  regions 

3)  effect  of  recombination  in  the  gas  (pressure  and  composition) 

4)  effects  of  the  energy  balance,  including  radiation  and  the  elevation 
of  the  electron  temperature 

The  list  can  be  rather  long.   We  closely  examined  the  energy  balance  because 
the  instabilities  are  believed  to  be  thermal  disturbances  (i.e.,  the  tempera- 
ture becomes  non-homogeneous)  and  recombination  in  the  gas  because  it  closely 
reflects  the  effect  of  pressure  and  composition.   For  simplicity,  a  two- 
dimensional  (infinite  parallel  plate)  geometry  was  considered,  see  Fig.  1. 

When  E   is  constant,  7* J  =  0   implies  that  J  =  J(y)   or  that 
n  =  n(y)   because  J  =  0  ;  moreover,   7  x  E  =  0  implies  that  E  =  E(y)  . 

y  y 

The  quasineutral  equations  for  steady  state  may  be  written  as   ' 
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Figure  1.   Geometry  of  Positive  Column 
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at     y  =   0  3n/3y  =   3T/3y  =   0 

at     y  =  +  a       n  =  0     and     T  =   T 
—  o 

In  the  above  formulation  x  is  the  axial  coordinate  and  y  the  interplate 
coordinate,  n  the  charge  concentration  in  the  ambipolar  region,   T   the 
electron  temperature  and  T  the  gas/ ion  temperature,  N  is  the  concentra- 
tion of  the  neutrals,   D   the  coefficient  of  ambipolar  diffusion,   K  the 

thermal  conductivity,   m   the  electron  mobility,   v.   and  a   the  ionization 

e  i 

and  recombination  coefficients  respectively,   E   the  axial  electric  field, 
v    the  collision  frequency  between  electrons  and  neutrals,   e  the  elec- 
tronic charge,   p   the  pressure,   and  m   and  m        the  mass  of  the  electrons 
and  neutral  respectively.   There  are  several  assumptions  inherent  in  this 
formulation,  namely  that  N>>n  ,  T  >T  ,  and  that  p  =  constant  .   In  addition, 

we  are  using  a  symmetry  condition  in  the  boundary  conditions  and  we  are 
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neglecting  of  wall  space  charges.   Glow  discharge  data   '   ,  in  addition, 

indicate  that   E   as  well  as   T   appear  to  be  constant  in  the  discharge  at 

least  before  constricting.   The  electron  energy  equation,  Eqn.  3,  should  have  to 

reflect  the  energy  dependance  of  the  collision  cross  sections  and,  in  addition, 


non-elastic  collisions  —  since  we  are  interested  in  laser  operation. 
However,  we  shall  not  deal  with  Eqns .  (3)  and  (4). 

Using  simple  kinetic  theory  to  describe  the  gas  temperature  dependance 
of  the  coefficients   (D   ,  p   ,  K  ,  etc.)   and  defining  a  stationary  state 
where 

v.n*  -  an*2  =  0  (5) 

i 

or       n*  =  v./a 

we  can  proceed  to  nondimensionalize  Eqns.  (1)  and  (2).   We  shall,  however, 

not  use  the  plate  spacing  to  render  y  dimensionless,  but  introduce  two 

characteristic  lengths,  L   and  L_   .   The  gas  temperature  is  made  non- 

n        T 

dimensional  by  the  use  of  a  wall  temperature  T   .   The  resulting  equations 


are: 


-  d/dy  (T1/2  0)  =  n  -  n2  (6) 


T   D 
V   xo  o 


d/dy  (T3/2  ||)  =  nT  (8) 


T  K 

V  en*u  ET 
eo  L 

It  seems  reasonable  to  expect  that  in  a  conventional  (unconstricted) 
glow,  the  two  characteristic  lengths  be  about  equal  to  each  other  and  that 

L  /a  x  0(1)   and  L  /a  «  0(1)  (10) 

n  T 


This  expectation  is  also  required  by  the  eigenvalue  nature  of  Eqns.  (6)  and 

(8).   It  may  be  shown  (see  Appendix  B) ,  that  except  for  the  equilibrium  limit, 

Eqn.  (6)  is  an  eigenvalue  equation,  (the  classical  solution  being  a  cosine 

when  no  recombination  is  present  and  when  T  «  constant) .   The  nature  of 

Eqn.  (6)  is  discussed  further  at  the  end  of  this  section. 

In  the  light  of  Eqn.  (10) ,  it  is  interesting  to  look  at  some  actual 

values  of  L   and  L_  .   Particularly  extensive  information  for  helium  is 
n        T 

reported  in  Ref.  16,  and  this  is  used  in  estimating  the  values  shown  in 
Fig.  2,  which  are  limited  to  the  range  of  applicability  of  the  equations 
in  Ref.  16.  An  appropriate  form  for  the  coefficients  is  obtained  by  intro- 
ducing the  gas  density  N   , 


T   (D   /N   ) 
Ln  No  =  V~2T     tv./N)  ™ 

v  O  1       o 


/  T     K     a 

LT  N      =    (E/N   )-1  W-, 2-pTT j— T  (12) 

To  o  y  e(y        N   ) (u . /N  ) 


The  above  forms  show  "similarity"  since  LN  is  a  function  of  E/N  only. 
Examining  Fig.  2,  it  is  evident  that  L  N  is  not  equal  to  L  N  ,  certainly 


no 


t  for  large  E/N  ,  wheareas  L^  follows  closely  data  reported  for  helium 
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in  the  literature   '    (the  fact  that  we  are  looking  at  parallel  plates  and 

the  data  is  for  tubes  is  immaterial) .   Since  in  a  discharge  tube  n  and   T 

share  a  common  diameter,  it  is  necessary  to  reconcile  the  discrepancy  of  the 

two  curves . 

14 
It  is  possible,  that  the  assumption   which  brings   T   into  Eqn.  (1) 

and,  therefore,  into  Eqn  (11)  is  invalid.   But,  even  if  we  were  to  eliminate 


..  o> 


..  00 


(N/3)  901 


-.  h- 


-*  t 


..  5f 


..  ro 


C\J 


the  factor  T  /T   from  the  equation,  the  curves  would  diverge  at  the  higher 
values  of  E/N  .   Constriction  of  the  positive  column  is  another  possibility 
but  in  this  case  L   would  be  expected  to  be  greater  than  L   for  any  given 
E/N  and  just  the  opposite  seems  to  be  true.   Uncertainties  in  the  estimate 
of  a  ,  the  recombination  coefficient,  would  affect  the  value  of  L   but 
this  estimate  is  the  one  that  follows  nicely  the  experimental  values. 

The  method  of  reconciliation  of  the  two  lengths  proposed  here  is  based 
on  the  alternative  that  the  equations  previously  written  are  incomplete  and 
that,  in  particular,   v.   in  Eqns .  (1)  and  (7)  should  be  replaced  by  a 
parameter   (C„)   which  makes  up  for  the  difference  in  the  two  curves.   This 
matter  is  discussed  in  Appendix  A  of  this  report. 

Consider  now  Eqn.  (1)  with  constant  coefficients.   We  proceed  here  to 
nondimensionalize  it  in  a  conventional  way 

-  D  n"  =  v.  n  -  an2  (13) 

a       i 

n'(o)  =  0  and  n(+  a)  =  0 

Let  n(o)  =  n   and  £  =  an  /v. 
o  o   i 


Hence  we  get 


-  a"  =  n  -  en2  d*) 

n'  (o)  =  0  and  n(o)  -  1 


n(+  a/v  D  /v.)  =  0 


Now  the  Shottky  solution  is  one  for  which  e  =  0  ,  namely, 


(15) 


and     n  =  n  cos  (tt  y/a)  (16) 

o     Z 


It  can  be  seen  therefore  that  the  characteristic  length  L  should  be 


L  - 


tt/2 


(17) 


On  the  other  hand,  for  the  equilibrium  limit  e  =  1  and  the  equation 

'2 


-  n  =  n  -  n 


(18) 


is  of  the  "Boundary  Layer"  type,  i.e.,  it  is  not  an  eigenvalue  problem,  see 
Ref .  14  pp  149-51  for  the  equivalent  problem  with  three-body  recombination. 
The  solution  to  Eqn.  18  in  an  implicit  form  is 


y  =  log 


|+m+|)1/: 


/!-rn+i)1/2 


•  -  1.317 


(19) 


or  explicitly 


n  =  —  n 
2   o 


1 

_ 

e 

(y 

+ 

i, 

.317)" 

1 

+ 

e 

(y 

+ 

l 

317) 

(20) 


which  is  indeed  of  the  boundary-layer  type.   The  intermediate  cases, 
0  <  e  <  1  ,  are  represented  by  elliptic  functions,  and  this  matter  is  dis- 
cussed in  Appendix  B  of  this  report. 
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IV.   CONCLUSIONS 

Some  deficiencies  of  the  Schottky  solution  have  been  re-examined  with 
the  intent  of  establishing  a  more  useful  solution  for  the  description  of 
the  transition  to  the  glow  collapse.   At  least  two  major,  recent  attempts  in 
the  literature  (Refs.  16  and  17)  base  the  steady  solution  for  the  charge 

profile  on  the  Schottky  result  and  these  attempts  have  not  succeeded  in  pre- 
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senting  a  complete  formulation  of  the  problem.   The  Technion  approach   is 

to  base  the  perturbation  on  a  physically  complete  solution  of  the  ambipolar 

charge  profile.   This  profile  must  include  the  ingredients  of  constriction 

and  of  the  elevation  of  the  electron  temperature. 

As  shown  in  Appendix  A,  the  characteristic  length  scales  for  ambipolar 

diffusion  and  heat  diffusion  (conduction  only)  may  be  reconciled  if   v.    in 

Eqn.  (7)  is  replaced  by  Cn  where  C„  >  v.    and   C_   is  a  function  of 
^  2         2    10        2 

E/N   (i.e.,  the  dependance  given  by  the  difference  of  the  two  curves  in  Fig. 
2).   What  appeared  to  be  a  bonus,  namely,  a  insight  into  the  mathematical 
description  of  striations,  needs  further  study  in  order  for  it  to  make  more 
physical  sense.   The  eigenvalue  nature  of  the  diffusion  equation  is  not 
charged  by  the  recombination  term  except  at  the  limit  when  z   =  1.0. 
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APPENDIX  A  -  PROBLEM  REFORMULATION 

Equation  (1)  in  the  text  presupposes  steady  flow  along  with  n  =  n(y)  . 
Let  us  relax  these  suppositions  and  go  back  to  the  species  flux  equations 
which  may  be  written  as 

a  D 

~  +  V-  [n(-  m  E-fVn)]-4  (Al) 

at  en 

a  D. 

f±+   V-  [n(u  I  -  -±  V  n)]  =  h  (A2) 

dt         in 

Combining  and  eliminating  the  electric  field  we  get  the  ambipolar  diffusion 
equation  in  the  usual  manner  which  for  constant  coefficients  is 

3n/3t  -  D  V2n  =  n  (A3) 

a 

or    3n/3t  -  D&  f-M*-  +  -Mh  -  v  n  -  an2  (A4) 

\3x    3y  / 


The  above  assumes  two-body  recombination  as  well.   Now,  using  Eqn.  (A3)  to 

2 

eliminate  V  n  from  Eqn.  (A2)  yields 

»     u  .  V •  nE 

3F  +  I-  D./D  =  n  (A5) 

1  a 

But  since  7*E  =  0   in  the  ambipolar  region 

H+1-"S/P  B-7n-°  (A6) 

i  a 

°r    3n/8C  +  (l  -"d./dJ  (EL  t  +  Ey  t) '  *  Vi  »  "  ™  (A7) 
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Let  us  work  with  Eqns.  (A4)  and  (A7)  in  linear  form,  i.e.,  neglecting 
recombination.   If 

n(x,y,t)  =  f(y)  g(x,t)  (A8) 

Then,  substituting  into  Eqn.  (A4) 

f  gt  -  Da  (fg^  +  gf")  =  v.fg  -  ...  (A9) 


d2f  32 

where         f"  =  — =■     g   =  — §■  and  g  =  3g/3t 

dy       **   3x        Z 
Dividing  by  fg  both  sides  of  (A9) 


1/8  («t  "  Da  ^   -  t   f"  =  Vi  (A10) 


At  a  given  E/N  ,  v.   is  constant  so  we  may  presume  that 

1/g  (g  -  D  g  )  =  constant  (or  -  C  )  (All) 

t      3   XX  X 


and  -  t2-  f "  =   constant    (or   CJ  (A12) 


with  -   C.    +  C0   =   v.  (A13) 

1         2  1 


The   constants      C,      and     C_      are  positive  and  defined   so   that      C_    >   v. 
12  2—i 

as  required  by  Eqn.  (A12) ,  which  now  replaces  Eqn.  (6)  since 

-  D  f"  •  f  C.  (A14) 

a         I 


and  Lr    =   7  D  /C. 

t      a  A 
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Equation  (A14)  is  the  new  Schottky  equation.   The  constant  C   appears  in 
the  companion  partial  differential  equation 

«t  ■  °a  gxx  +  Cl  «  =  °  (A15) 

Now  look  at  Eq.  (A7) .   At  the  centerline   3n/8y  =  0  and  f(o)  =  f 

o 

so  that 

acy  =  0  at*  (l  -  d1/dJ-el  te  '  vi  n  '  •••  (A16) 

Introducing  n  =  fg  we  obtain 

«t  fo  +  (l  -  D./dJ  EL  *x  £o  "  Vi  £0«  CA17) 

We  can  cancel  f   in  the  linear  equation  above  and  obtain  an  equation 
in  which  no  function  of  y  appears 

»t  +   (l  -\/dJ  EL  «,  ■  »1«  (A18) 


At  this  point  it  is  worthwhile  to  eliminate  g   between  Eqns.  (A18)  and 
(A15)  and  obtain 


Da  8xx  "  Cl  g  =  V  "   1  -  D./D   V 


Ui 


x 

\     i"  a/ 


"A 


Da  ^  +  (1  -  D./D  )  ^  ■  ('Ji  +  C1)S 


i  a 
°r       Sxx   \D  -  V±J   Sx  ""  Da  g  (A19) 
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Equation  (A19)  can  be  solved  with 


g  =  T(t)  e  aX  (A20) 


where  a  is  found  from 


-  C  +  \li2  +   4C9/D  y.E_ 

a    i 


Note  that   D.  <  D   and  that  c,      is  anyway  a  positive  parameter. 
Equation  (A21)  indicates  that  the  roots  are  real  and  there  is  no  periodic 
behavior  in  x  —  that  striations  are  not  to  be  found  with  this  description. 

In  order  for  n   to  oscillate  in  x  ,  we  must  have  a  situation  in  which 
new  constants  C  *  and  C~*  are  introduced  such  that 

C  *  -  C2*  =  v±  (A22) 

§t  "  °a  «»  =  °1*  S  (A23) 

and     -  D  f"  =  -  C  *  f  (A24) 

a        z 

Here,  Eqn.  (A24)  will  result  in  (non-periodic)  hyperbolic  functions  in  y  and 
Eqn.  (A24)  in  oscillatory  functions  in  x  provided  that 

4C2*    h±\  V 

ir-  '^T-\  (A25) 

As  before  we  would  like  to  have 

Cl*  >  Vi   (C2*  =  Cl*  "  Vi}  (A26) 

so  as  to  match  L   to  L_  . 
n       I 
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The  solution  to  Eqn.  (A24)  does  not  represent  the  usual  diffusion 
controlled  situation  because  in  order  to  meet  the  symmetry  condition  at  the 
center line  we  must  give  up  the  physically  sensible  boundary  conditions  at  the 
wall.   Perhaps  variable  coefficients  and  the  inclusion  of  recombination  can 
remedy  this  problem  with  the  boundary  conditions.   Or  perhaps  we  should  be 
willing  to  give  up  axial  symmetry  and  consider  some  sort  of  contraction 
asymmetry  which  would  likely  have  to  include  some  reconbination  because  of 
the  higher  pressures.   This,  however  shall  not  be  pursued  in  this  report. 
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APPENDIX  B  -  INFLUENCE  OF  TWO-BODY  RECOMBINATION  OF  SCHOTTKY  SOLUTION 

The  diffusion  equation  given  in  the  main  text  (Eqn.  14)  is  not  suitable 
for  elliptic  funciton  representation,  therefore,  a  different  form  is  developed 
below. 

The  diffusion  equation  with  constant  coefficients  is  (without  considera- 
tion of  whether  C?   shall  be  used) 

a2  ? 

-  D  2-~  =  v .  n  -  an  (Bl) 

Sl    .    2.  i 

dy 

2 
Defining  a  stationary  state  v.n*  -  an*  =  0  ,  or 

n*  =  v./a  (b2) 

We  can  nondimensionalize  the  equation  with 

n  =  n/n*  and  y  =  y/L 


n 


Obtaining 

_  <L_|  =  £  _  ~2  (B3) 

dy 

n'(o)  =  0  and  n  =  6  ,   0  <_   8  <_  1 

n(a/L  )  =  0 
n 

It  is  clear  that  the  boundary  condition  at  the  centerline  make  up  for  the 
disappearance  of  e      in  Eqn.  (B3) .   Now  Eqn.  (B3)  can  be  integrated  once  to 
yield 


to  .^(2  -3  .  .2)  _(2g3_B* 
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\jf  f  d"  f  <U  (B5) 

V   ft7   Vn3  _  2  n2  _  s3  +  |  62   ^ 

The  left-hand- side  of  the  above  equation  can  be  represented  in  terms  of  ellip- 
tic functions.   Let 

(n3  -  |  n2  -  32  +  |  e2)  =  (n  -  a)  (n  -  b)  (n  -  c)  (B6) 

therefore 


a 


}-»)+#-•)$♦»' 


b  =  3         (0  <  3  <  1)  >     (B7) 


j-»)-V(J- »)(!  +  * 

C  =   r 


where  a  >  b  _>  n  _>  c 

With  information  on  a,  b,  and  c  we  can  enter  any  table  for  elliptic 
functions  (see  for  example  Ref.  18  or  19).   Using  Byrd  and  Friedman  notation, 
formula  234.00,  we  have 


b 

/ 


V(a  -  n)(b  -  n)(n  -  c) 
h 

where  2         2   b  -  c 

8    V*  ~  c  '       '  a  -  c 


=  gF  (*>,k)  (B8) 
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and  „  -  sin"1  J<*  -  ^  Z   ?> 

\r(b  -  c)  (a  -  n) 

F(^,k)   is  tabulated.   Note  that  for  6  =  1  ,  we  get  a  double  root   (a  =  b) 
and  the  equation  is  not  an  elliptic  integral  (though  its  solution  can  still  be 
found  from  the  Tables) .   When  8=1  we  get  the  boundary-layer  behavior 
described  by  Eqn.  (19)  or  (20)  in  the  main  text. 

Except  for   6=1,  Eqn.  B3  represents  an  eigenvalue  problem.   In  order 
to  check  this,  we  solve  for  the  value  of  the  integral  at  the  wall  where 
n  =  0  .   Thus 


.  -1   /(a  -  c)b 
'  =  Sin    '(b  -  c)a 


For  very  small  values  of   6  ,  we  obtain  the  cosine  solution,  namely, 

X   =  (^/2)2    where    L  =  a/X1/2 

n 


As   6  increases,  we  obtain  progressively  larger  values  until  at   6=1  we 
get  A  =  °°  )  .   Figure  Bl  shows  this  behavior  for  the  eigenvalue.   The  results 
shown  in  Fig.  Bl  proved  identical  to  results  from  a  computer  program  at  the 
Technion  (for  the  combined  diffusion  and  heat  dissipation) . 

Increasingly  higher  degrees  of  recombination  flatten  out  the  density 
profile  at  the  center  and  increase  its  slope  at  the  walls.   The  radial  elec- 
tric field  is  thus  enhanced.   The  variable  coefficient  solution  is  not  antici- 
pated to  change  the  above  results  in  any  major  way. 
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Figure  Bl.   Eigenvalue  for  Increasing  Degrees  of  Recombination 
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